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Abstract. We present experiments on the luminescence of excitons confined in 
a potential trap at milli-Kelvin bath temperatures under cw-excitation. They 
reveal several distinct features like a kink in the dependence of the total integrated 
luminescence intensity on excitation laser power and a bimodal distribution of the 
spatially resolved luminescence. Furthermore, we discuss the present state of the 
theoretical description of Bose-Einstein condensation of excitons with respect to 
signatures of a condensate in the luminescence. The comparison of the experimental 
data with theoretical results with respect to the spatially resolved as well as the 
integrated luminescence intensity shows the necessity of taking into account a Bose- 
Einstein condensed excitonic phase in order to understand the behaviour of the trapped 
excitons. 
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1. Introduction 

Almost 50 years ago, excitons [?, ?] have been suggested as particularly interesting 
candidates for Bose-Einstein condensation (BEC), as they consist of an electron and a 
hole in a semiconductor, both fermions bound to form a bosonic excitation and thus 
resembling most closely neutral atoms of usual matter. Due to their rather small mass 
comparable to the free electron mass, it was speculated that for exciton densities of the 
order of 10 18 cm" 3 - easily achievable by absorption of photons - critical temperatures 
of some 10 K may be reached. 

Due to their unique properties, the excitons of the so-called yellow series in the 
semiconductor cuprous oxide (CU2O) are still considered the most promising candidates 
for excitonic BEC [?, ?, ?]. This is related to the large binding energy of 150 meV, which 
shifts the Mott density to 3 • 10 18 cm -3 at cryogenic temperatures [?, ?]. Made up from 
doubly degenerate valence and conduction bands, the ground state of this series splits 
into the triply degenerate orthoexciton and the nondegenerate paraexciton, which is the 
energetically lowest exciton state, lying A = 12.1 meV below the orthoexciton states. 
Due to the positive parity of the bands, the orthoexciton is only weakly optically allowed 
(quadrupole transition with oscillator strength 3 ■ 10" 9 [?]), while the paraexciton as a 
pure triplet state with respect to the electron and hole spins [?] is optically forbidden 
in all orders. Its decay is only possible via an odd parity optical phonon resulting in a 
long lifetime in the microseconds range during which thermodynamic quasi-equilibrium 
may be reached. 

As in all physical systems for which BEC has been demonstrated up to now, excitons 
should be confined in a potential trap. This has the advantage that (i) the diffusion 
process, which reduces the exciton density, is suppressed and (ii) the critical number of 
particles required for the phase transition decreases much faster with temperature than 
in free space. The critical particle number is given by 



where VLq is the average oscillator frequency of the trapping potential and ( denotes the 
Riemann Zeta function [?] . 

Despite the promising properties mentioned above, all previous studies to create 
a dense gas of excitons in CU2O either in a bulk crystal or in a potential trap did not 
demonstrate conclusively excitonic BEC [?, ?, ?, ?, ?, ?]. The main reason for this 
failure turned out to be the existence of a very efficient exciton-exciton annihilation 
process that sets in at high exciton densities and whose rate scales with the square of 
the exciton density n 

dn 9 

*=-<•" • M 

Due to the rather large value of a of the order of 10~ 16 cm 3 ns _1 , this process was believed 
to prevent the formation of a BEC of paraexcitons. 




(1) 
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However, in recent experiments using pulsed excitation, we have found in agreement 
with earlier studies [?] that this process can be damped in shallow potential traps by 
almost two orders of magnitude [?]. While this allowed the accumulation of large exciton 
numbers, the temperatures in these experiments were still too high to undercut the 
critical temperature for Bose-Einstein condensation, which was predicted to require a 
bath temperature as low as 100 mK at the experimental conditions. Furthermore, the 
high pulse energies needed to create sufficient exciton numbers resulted in a substantial 
heating of the exciton gas during the laser pulse. Our strategy to overcome these 
problems was twofold: On the one hand, we reduced the temperature of the He bath 
as much as possible by preventing thermal radiation from the surrounding to reach the 
sample. In this way we reached a minimum temperature of 35 mK for zero incident laser 
power, which increases at laser powers of about 1 mW to 250 mK within one hour. On 
the other hand, we switched over to cw-excitation, for which the possibility of creating 
large paraexciton numbers has been demonstrated recently [?]. By reducing the thermal 
load, we were able to get down to effective exciton temperatures as low as 200 mK at 
low excitation power. The main advantage of cw-excitation, however, is the possibility 
to achieve a quasi-equilibrium situation in which the decay rate of excitons is cancelled 
by the formation and relaxation rate of the species in the trap. This allows to drive the 
system through the possibly existing various phases by simply changing the excitation 
power of the driving laser, as was exemplified by the condensation experiments of exciton 
polaritons in semiconductor microcavities (for a review, see [?]). In contrast to polariton 
systems, the lifetime broadening for paraexcitons is extremely small due to their long 
lifetime of the order of [\s. Therefore, we can neglect any damping effects due to the 
excitonic decay in the spectra [?, ?]. 

The paper is organised as follows: In section [2j we briefly sketch the experimental 
setup and present a typical set of experimental results which reflect a variety of 
experimental conditions obtained by changing the bath temperature and the excitation 
conditions. We concentrate in particular on the spatial profiles of the luminescence 
intensity and on the totally integrated intensity in dependence on the excitation 
power. These results are analysed theoretically in section [3] calculating the excitonic 
luminescence of an interacting Bose gas on a mean field level. Thereby, we start with 
the usual assumption of global equilibrium and then extend the theory to the case of 
excitons in local thermodynamic equilibrium. Finally, we compare the theoretical results 
with the measurements and show that there is excellent qualitative agreement if we take 
into account the occurrence of a Bose-Einstein condensate of excitons. Section [4] gives 
our conclusion and an outlook to further experiments. 
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2. Experiment 

2.1. Experimental setup 

For the studies at subkelvin temperatures, we used the same experimental setup as 
reported previously [?], but implemented a narrow-band tunable dye laser (Coherent 
CR599, laser dye Rhodamin 6G) pumped by a 5W green solid state laser (Verdi 5), see 
figure [TJ The laser power was stabilised by a closed feedback loop to within 1%, the laser 
frequency and line width (< 0.5 GHz) were measured with a wavemeter (High Finesse 
WS7, resolution 60 MHz). In order to enhance both the spectral and spatial resolution, 
we employed a fourfold magnification optical imaging system between the spectrometer 
exit slit and the detector. 

For the experiments we used natural cuprous oxide crystals originally found in 
Namibia in the form of millimetre sized cubic specimens with well defined facets (see 
[?] for details). The quality of these samples was checked according to a low defect 
density, leading to long paraexciton lifetimes up to 1 \xs. For such samples, previous 
high resolution absorption measurements in a magnetic field revealed a paraexciton line 
width as narrow as 80neV, demonstrating their extremely high quality [?]. 

The potential trap for the confinement of the exciton gas was made by the well- 
known Hertzian stress technique [?, ?, ?, ?], where a spherical stressor made of glass 
(radius 7.75 mm) is pressed with a force F against a flat surface of the crystal along a 
direction which we denote as z-direction (figure [I]). As a result, a confining potential is 
generated in which the energies of ortho- and paraexcitons are lowered compared to the 
bulk. This potential V eX t(x, y, z) can be calculated from the known strain parameters 
of the yellow exciton states (for a recent calculation with refined parameters see [?]). 
To achieve agreement between the calculated potential profiles and the experimentally 
measured low energy border lines of the spatio-spectral images, we had to use as stressor 
radius a value 50 % larger than the nominal one [?]. For a simple description, we 
decompose the potential trap and give in table [T] the parameters of the harmonic 
oscillator potential along z, V(z) = an (z— Zq) 2 — Vq, and of the two-dimensional harmonic 
oscillator normal to z in the x,y— plane, V(x,y) = a±(x 2 + y 2 ). 

As resonant creation process of the excitons we used the indirect absorption process 
involving an odd parity optical -phonon into the orthoexcitons in the trap (see inset 
in figure [2]). These orthoexcitons quickly transform into paraexcitons. The laser beam 
was positioned about 100 |xm away from the trap centre in positive z-direction, i.e., 
away from the stressor lens. In order to confine the primarily created orthoexcitons in 
the trap and to avoid any excitation outside the trap, which may lead to losses of the 
excitons, we tuned the energy of the laser photons slightly (« 0.5 meV) below the onset 
of the phonon sideband in the bulk (at 2048.56 meV). By this process, about half of the 
incoming photons are transformed into excitons. 

To detect the primarily excited orthoexcitons for the calibration of the number of 
excited excitons, we have chosen the k- vector of the excitation beam along the direction 
of observation. 
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parameter 



value 



potential curvature ay 
potential curvature a± 
lifetime of paraexcitons 
ortho-para conversion rate 
two-body decay rate App 
two-body decay rate Aqq 
relaxation rate r re i 



0.1334 [leY/imi 2 
0.0733 [leV/iim 2 
650 ns 
0.2 ns" 1 f 

2-10" 18 cm 3 ns" 1 { 
4.9 -10 -17 cm 3 ns" 1 
6-10 7 s" 1 



potential trap minimum Vq 1.35 meV 



Table 1. Parameters for trapping potential and rate model, 
f taken from reference [?] 

| due to the lower strain, the two-body decay rate App is reduced by a factor of two 
compared to [?] 



In the experimental setup (figure [I]), the emitted light was imaged onto the entrance 
slit of a high resolution triple spectrograph (T64000, Jobin Yvon) usable either in 
subtractive or additive mode with a diffraction limited spatial resolution of the order 
of 10 [im. The astigmatism of the spectrograph was corrected by a cylindrical lens 
(focal length F = —1000 mm/ F = 1000 mm for subtractive and additive dispersion, 
respectively) in front of the entrance slit [?]. To obtain a ^-resolved spectrum I(z,u), 
along the direction of the applied strain, the luminescence from a small stripe of width 
2Ay centred in the trap was integrated along the x-direction perpendicular to z (see 
figure [T|. Detection was done either by an intensified CCD camera (Andor iStar) 
which could be gated with a minimum temporal resolution of 5 ns or with a nitrogen 
cooled CCD camera with high quantum efficiency (Andor Newton), which allowed long 
integration times. 

2.2. Application of the rate model to cw- excitation 

With some adaptions, the rate model developed in [?] can also be applied to describe 
the cw experiments. These concern the following: 

• Even under resonant excitation of excitons, a considerable density of unbound 
electron-hole pairs is generated [?]. This implies that the formation of an exciton 
from the hot electron-hole pair generated in the two-body decay (Auger-like process) 
takes a finite time comparable with other exciton relaxation times, e.g., the one 
related to phonon scattering. In a fast initial relaxation stage with a duration of 
several picoseconds [?], the carriers thermalize by LO and LA phonon scattering. 
For a strained crystal with an exciton potential trap, these electron-hole pairs will 
then undergo different relaxation scenarios depending on whether the strain causes 
a trapping potential also for the unbound electron- hole pairs or not. In the latter 
case, the electrons and holes will diffuse from the place of generation (the potential 
trap) into the whole crystal and form excitons which themselves will drift again 
into the trap. In the former case, the electron-hole pairs will stay inside their trap 
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Figure 1. Main features of the experimental setup. Excitons are created by a dye 
laser propagating along the ^-direction normal to the strain axis (z-direction). The 
laser was focused either directly into the trap or onto a spot close to the trap. In the 
latter case the excitons diffuse towards the trap due to the gradient force initiated by 
the potential trap. The emission out of the trap is monitored spatially resolved along 
the y- and z-direction, while integrating along the x-direction. 



and form a stable electron-hole plasma cloud with a density which is three orders 
of magnitude smaller than the exciton concentration according to the results in 
[?]. As derived in Appendix A, strain leads to a trapping potential for unbound 
electron-hole pairs which is similar to that for the paraexcitons. To take this effect 
into account, we included in the rate model the unbound electron-hole pairs with a 
total number of iV eh p. They are generated via the two-body decay of the ortho- and 
paraexcitons and recombine with a rate T rc N^ hp . The recombination rate should 
depend on temperature in the same way as the two-body decay of the excitons 
(ocT- 3 / 2 ). 

Due to heating and incomplete relaxation, the excitons may not cool down to 
bath temperature. We can describe this effect by assuming an effective exciton 
temperature T eS , a cooling time of 200 ns [?] and heating processes due to the 
energy release by ortho-para conversion Co = 10 meV per exciton, Auger-like 
two-body decay with Cxx = 2 eV per exciton pair, and a non-radiative decay of 
paraexcitons (guessed release Cp = 0.25 eV per exciton). The temperature rise 
due to this heating is characterized by a constant Ch ea t which was adjusted to 

C h eat = 7-10 9 K/J. 
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We thus obtain the following set of rate equations for the numbers of primarily 
excited orthoexcitons (O exc ), trapped ortho- and paraexcitons (Ot and P^, respectively), 
and unbound electron-hole pairs iVehp, and the effective exciton temperature T e g: 

d ^ L =N - G(t) - T rel ■ O exc - r OP ■ O cxc - T ■ O exc , (3) 



- Trci ■ Ocxc — Top • Ot — 2Aqo • 0% 



dt 



dPj 

~dT 



Aop . 0t .P t -T -O t + ^T rc Ni p , (4) 



Top ■ Ot + Top ■ O cxc — 2App ■ P^ — Tp ■ P^ 

-A p-o T -p T + ^r rc iv e 2 hp , (5) 

— - App ■ P| + A 00 ■ 2 T + A p -O t -P t - T rc iV c 2 hp - r ehp iV ehp , (6) 




0.01 



0.1 



1 10 100 

P L ( M W) 



1x10^ 



1x10 



Figure 2. Numbers of para- and orthoexcitons and unbound electron- hole pairs and 
effective exciton temperature as a function of excitation laser power as expected from 
the rate model. The inset schematically shows the laser excitation process to create 
the orthoexcitons, their conversion into paraexcitons and the luminescence processes 
to detect the excitons. 
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dT cS 



dt 



a 



heat 



Cxx (^4pp ■ Py + ^oo " 0\ + A Q p ■ Ot ■ Pt 



-CqFop{Ot + O exc ) + CpTp ■ P T ] — T coo i(T cS — T fin ) . (7) 



Here, Tfi n denotes the final temperature to which the exciton gas would relax without 
additional heating. In an equilibrium situation, this would correspond to the bath 
temperature. 

To simulate the cw-excitation, we solved the system of rate equations by assuming 
for G(t) a rectangular shaped excitation pulse of width AT L = 25 [xs and unit pulse 
area. The number of initially excited excitons iVo is given by 

A Pt 

N = ^AT L , (8) 

with the conversion factor from incident laser photons into primarily excited excitons 
A = 0.45. All other parameters are taken from reference [?]. Since we used the same 
sample under similar conditions, the calculation allows to obtain the number of para- and 
orthoexcitons in the trap in dependence on laser power, whereby we guess the accuracy 
to ±50%. In figure[2]we show a set of typical results. The final temperature T fin = 0.35 K 
was chosen to reproduce the experimentally observed temperature dependence in figure 
[5p. The results at high power levels do change somewhat due to the temperature 
dependence of the Auger process, but not more than 20% by varying T^ n from 0.2 
K to 0.5 K. 

The number of paraexcitons in the trap at low laser powers turns out to be 
determined by two parameters, the paraexciton lifetime 1/Tp and the fraction of 
absorbed photons A. From figure [2] one can estimate the number of paraexcitons in the 
trap to be about 8 • 10 6 for a laser power of 10 |J.W. The number of unbound electron- 
hole pairs would be around 3 ■ 10 5 . At power levels of lmW, which was the maximum 
laser power used in the mK experiments, paraexciton numbers of 2 ■ 10 8 can be realised 
and the number of electron-hole pairs increases to 10 7 . Assuming thermal equilibrium, 

ehp 



this would correspond to densities of about n-p = 2 ■ 10 16 cm 3 and n c \ ip = 10 14 cm 3 . 



2.3. Experimental results 

2.3.1. Low excitation power The first experimental results we want to show were 
obtained at rather low laser excitation powers in the range between 1 and 5 |J.W with 
high spectral resolution (figure [3]). Excitation was performed via the phonon-assisted 
absorption of the orthoexcitons with energy slightly below the bulk exciton gap as 



described in section 2.1 Orthoexcitons were converted rapidly into paraexcitons and 
relax down to the bottom of the trap. 

From the z-resolved spectra (panels A) we obtained the ^-integrated spectra (panels 
B) and the z-profiles (panels C). As already could be seen in the resolved spectra, there 
is no change in both the spectral and spatial distribution, only an increase in the overall 
intensity which is almost proportional to the laser power. This demonstrates that we are 
in the linear excitation regime and that bimolecular decay processes are not important. 
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Figure 3. Spectra and spatial luminescence profiles at low excitation power Pl with 
Pl = 1-5 M-W (upper row), Pl = 3|xW (middle row), and Pl = 5|J.W (lower row). T s 
denotes the temperature extracted from the high-energy tail of the spectra. 



The high-energy tail of the spectra can be fitted quite well by a Boltzmann distribution 
with an effective exciton temperature of about 0.4 K. A two-dimensional spatial image 
of the exciton cloud was obtained by using the spectrometer in the subtractive mode, the 
output stage at zero wavelength and by setting the intermediate slit of the subtractive 
stage to just let pass the emission from the paraexcitons in the trap. A typical example 
is shown in figure |4| The spatial profiles are described by simple Gaussian distributions 

I(y, z) = n exp (-v 2 /**) exp (- z 7^) (9) 



with width parameters a± and cry of about 20 llhi. 

According to the theoretical predictions in thermal equilibrium (see section [3]), the 
width parameters are related to the curvature of the trap a and temperature T by 



OX 



(10) 



However, the ratio of the resulting width parameters does not fit the ratio of the potential 
curvatures as obtained from the strain calculation. It comes out about 30% larger. As 
origin of this broadening we identified a small vibrational motion of the sample in y- 
direction of the order of 10 (xm which could be directly observed by imaging a submicron 
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sized needle. For a quantitative theoretical model of the light emission of the excitons, 
this effect has to be taken into account. 

Using the curvature parallel to the z-direction from table 1, and subtracting the 
spatial blurring due to diffraction of our optical imaging setup (7p,m), we obtain from 
the spatial profiles temperatures of about 0.45 K, in good agreement with those from 
the spectra. On a first sight, we would conclude that the excitons are in a state of 
quasi-thermal equilibrium, albeit with a temperature which is more than an order of 
magnitude larger than the bath temperature, which was in these experiments around 
38 niK. There are two possible explanations of this difference: 

(i) The heating of the exciton cloud due to the excess kinetic energy of about 8 meV 
in the ortho-para conversion process or due to the exciton-exciton annihilation 
process. The excess electronic energy is in both cases converted into phonons and 
finally increases the temperature of the lattice, which is unknown in the present 
experiment, however. This heat has to be transferred into the surrounding He bath, 
where it would lead to an increase in the temperature of the mixing chamber, which 
could be measured during the experiment. 

Indeed we observe this heating effect at higher laser powers. As shown in figure 
[5j the effective exciton temperature and the temperature of the mixing chamber 
rise from their initial values very steeply above P L = 1000 p.W. Presuming a well- 
defined relation between the lattice temperature of the crystal and that of the 
mixing chamber, we have to conclude from these data that, at low laser powers, 
the heating effect is not present and the lattice temperature should be close to that 
of the He-bath at very low laser powers. This conclusion is in agreement with the 
results of the rate model presented in the foregoing section, too. There the heating 
of the exciton cloud is starting at laser powers well above 100 p.W and is not present 
at low powers. 

(ii) The discrepancy of the lattice temperature Ti att and that of the exciton cloud may 




y (P' xe ') spatial position (pm) 



Figure 4. left panel: Spatial image of the exciton cloud at low excitation power of 
Pl = 30 (J.W. Right panel: Marginal intensity distributions along z (blue circles) and 
y (red circles) with fits by gaussians according to equation [9] 
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result because the excitons do not reach global thermodynamic equilibrium during 
their finite lifetime. This standard point of view means that, due to elastic exciton- 
exciton scattering, a quasi-thermal equilibrium state is established very fast with 
a temperature above that of the lattice. Hereafter, the exciton gas cools down to 
the lattice temperature by emission of acoustic phonons. The time scale of this 
process increases with lower lattice temperature. In an external potential provided 
by the trap, the excitons are additionally driven by force and drift terms into 
the potential minimum, as described by the Boltzmann equation [?]. The detailed 
balance of these processes is determined by the interplay between the various exciton 
relaxation mechanisms. 

In a recent paper, we published a detailed numerical simulation by solving the 
Boltzmann transport equation for excitons in a potential trap under pulsed optical 
excitation considering phonon relaxation and exciton-exciton annihilation [?] . The 
simulation indeed showed a strong non-equilibrium situation. While the local 
exciton energy distribution reached, within the exciton lifetime of several 100 ns, 
an equilibrium with the lattice down to temperatures above 300 mK, the spatial 
distribution, which is governed by the force and drift terms, remained much broader 
than in thermal equilibrium. For the case of still lower lattice temperature, which 
was not considered in [?], we obtain locally thermodynamic quasi-equilibrium also 
at temperatures down to 50 mK by including elastic exciton-exciton scattering as 
is shown in detail in Appendix B Albeit in this range the exciton temperature does 
not come down to the lattice temperature within the exciton life time, it is much 
lower than the observed values of T s . According to these results, in our experiments 
we should expect a non-equilibrium situation in which the excitons locally are in 




photon energy (pixel) 



Figure 5. Panel A: Series of spectra taken at increasing laser power Pl from 5 (xW to 
3 mW in a logarithmic intensity scale. The spectral position is given in CCD-pixels, 
1 pixel corresponds to 2.5 [ie V. Panel B: Spectral temperature T s as obtained from a 
Boltzmann fit of the high-energy tail of the spectra shown in panel A (left ordinate) 
and the temperature of the mixing chamber Tu e (right ordinate) . 
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quasi-thermal equilibrium with a temperature Tx, which we assume to be the same 
at all points in the trap. Globally, the exciton distribution will be quite different, 
but still can be described approximately by a Gaussian dependence on position. 

Considering the effect of such a situation on the experimentally measured 
luminescence spectra, we recall that only excitons near k participate in the light 
emission, where ko is the wave vector of the intersection of the photon and exciton 
dispersions. Its modulus is given by ko = E g xn/hc, with E gX being the excitonic 
band gap, n the refraction index, and c the vacuum velocity of light. This means, we 
observe the exciton distribution function / eq (k, r) only at one point in local A;-space. 
Furthermore, due to energy conservation, at a certain spectral position Tluj we only 
observe those excitons, which are at that spatial position in the trap, where the trap 
potential V ext (r) = %uj — E — e . Here, £"0 is the energy of the trap bottom and e is 
the kinetic energy of the excitons at ko- If we have a spatial distribution 



which suggests a Boltzmann distribution with an effective temperature T s = aa 2 /k-Q- 
Then the spectral distribution of luminescence intensity in reality describes the spatial 
distribution and vice versa. Therefore, it has to result in the same effective exciton 
temperature, which is also a consistency check for the trap potential. Consequently, 
we will designate this temperature as the spatial temperature T s of the exciton cloud. 
Obviously, T s corresponds to T eff of the rate model, cf. equation ([7]). On the other hand, 
it is not possible to extract the local effective exciton temperature Tx from the intensity 
of the zero-phonon transition. This would require a spectroscopic technique which is 
sensitive to the local exciton distribution like the transitions between IS and 2P exciton 
states (see [?, ?]). Hence, at the present stage of measurements, all we can state is that 
Tx may differ considerably from both the lattice temperature, which will set a minimum 
value, and the spatial temperature, which will set a maximum value, as in the latter 
case we would have global equilibrium. 

2.3.2. Power dependent experiments In the following we will present a series of 
measurements under controlled conditions of excitation and cooling times such that 
we are able to influence the thermodynamic state of the system in order to optimise 
any possible phase transition. The laser power ranged from below 1 \xSN to about 
700 [iW. For each measurement we noted the temperature of the mixing chamber (He- 
bath temperature T bath ) at the beginning and at the end. Between the measurements, 
the laser beam was blocked, so that the crystal could cool down (see figure |6|. 

In order to obtain the best signal-to-noise ratio, the ICCD camera was operated 
in the photon counting mode. Here, the CCD-chip is read out after a short exposure 
time (T) and the positions of the single photon peaks are determined and stored. A 





this yields a spectral distribution of the form 




(12) 
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single photon is identified, if the signal falls between a lower and an upper discriminator 

level, which were adjusted so that only single photon signals fall in this range. For a 

number iV acc of accumulations, the data are summed up in the memory. Advantages of 

this method are that signals due to cosmic rays are easily detected and eliminated and 

that one can estimate the error of the intensity signal due to the Poissonian statistics as 

the square root of the number of photons collected at each pixel. Due to the Poissonian 

statistics, the probability of detecting n photons is given by 
yn 

P n = — exp(-A), (13) 
n\ 

where A is the average photon number for each exposure. As we detect only single 
photons and throw away all higher photon number events, A should be much smaller 
than one. However, for moderate values, we can obtain the true number of detected 
photons, AiV acc from the measured number iVdet by applying the correction formula 

\N acc = -\n(l-^p) . (14) 

V -''ace' 

For the measurements of the z-resolved luminescence spectra shown in figure [7] we 
have Ti = 0.5 seconds and N acc = 400. In order to keep the maximum number of 
detected photons well below iV acc , we inserted neutral density filters before the entrance 
slit of the monochromator. Their transmission factors have been determined by direct 
intensity calibration as Dl = 0.314 and D2 = 0.0333. 

In figure [7j each measurement is numbered according to its position in time given in 
figure |6j For each image, we specify the laser power, the maximum number of counted 
photons, which allows to judge the signal-to-noise ratio, and the totally integrated 

95 i 1 1 1 1 1 




time (s) 

Figure 6. Time profiles of laser power (blue bars) and temperature of the 
mixing chamber (red diamonds), and spatial temperature (black diamonds) for the 
measurements shown in figure [7] The red line is a guide for the eye. 
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Figure 7. Series of z-resolved luminescence spectra for different excitation powers Pl 
and bath temperatures (see figure [6]). C max denotes the peak value and I t ot the totally 
integrated intensity. The numbers in parentheses denote the measurement numbers, 
cf. figure [6] 
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intensity. At first sight, all the images look very much the same, despite the three 
orders of magnitude increase in laser power. In particular, we observe neither the 
sharp spatial and spectral peak at the bottom of the trap, which is expected as a 
signature for BEC in case of an ideal Bose gas, nor the predicted flattening of the 
emission along the z-direction in case of an interacting Bose gas [?]. In this respect, the 
measurements shown are similar to those reported in all earlier work [?, ?]. A closer 
look, however, reveals subtle changes of the shape of the spectra depending on the laser 
power. For power levels below 30 |a.W, the ^-resolved spectrum of the exciton cloud 
shows no significant changes - here note especially measurements 9 and 21, which were 
taken under very different He-bath temperatures (see figure [6]). For power levels above 
30 |xW, the images change gradually by showing a spectral narrowing and a slight blue 
shift. These changes suddenly come to a stop between measurements 17 and 18, with the 
images from now on both spatially and spectrally much broader and resembling those 
taken at low powers. The latter data have been taken at almost the same laser power 
of 300 (J.W, but the He-bath temperature was significantly different. Measurement 17 
started at Tb a th = 42 mK. During the measurement, the laser power resulted in a heating 
up to Tbath — 49 mK. In contrast, measurement 18 already started at Tb at h = 45 mK 
and stopped at Tb a th = 50 mK. This temperature variation of the mixing chamber 
should directly reflect the actual lattice temperature T latt of the Cu 2 crystal. At the 
beginning of measurement 17, Ti att must have been quite close to that of the bath 
while at the start of measurement 18, T latt has only cooled down to midway between 
start and end of measurement 17. This difference in lattice temperature is reflected in a 
concomitant difference in the spectral temperature which is T s = 0.62 K for measurement 
17 and T s = 0.72 K for measurement 18. The totally integrated intensities of both 
measurements are expected to be similar due to the same laser power, however, that of 
measurement 18 is about 30% larger. 

We can resolve this puzzle by assuming that the difference between both 
measurements comes from the existence of a Bose-Einstein condensate of excitons 
in measurement 17, while measurement 18 describes a normal exciton cloud at 
temperatures above the critical temperatures for BEC for the actual number of excitons 
in the trap, which is the same in both measurements. Then the less intensity in 17 
would be the result of the decreased number of thermalised excitons, since according to 
[?] and anticipating the discussion in section |3j excitons in the condensate are in the 
ground state (k = 0), and, therefore, could not emit light with wave vector fc at all (as 
in a homogeneous system), or strongly suppressed (as would be the case in the trap). 

Thus, we propose the following scenario for measurement 17: In the beginning, we 
have a system with a large fraction of the excitons in a condensate, with a reduced 
luminescence intensity. As time increases, the crystal heats up and the condensate 
fraction is reduced. When the lattice temperature reaches the critical temperature, we 
have a situation just like in measurement 18. Therefore, after being integrated over the 
exposure time Ti, the ^-resolved spectra should look the same. In a first approximation, 
the image of 17 is a superposition of that of measurement 18 reduced by a factor 7 and 
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that of a system with a condensate. To demonstrate that this is indeed the case, we 
compare in figure [8] in the left panel the z-profiles for both measurements. By assuming 
7 = 0.6, the wings of both profiles do coincide exactly, which is clearly visible in the 
difference profile. The difference around the trap centre, however, has a characteristic, 
non-Gaussian line shape resembling that of a condensate (compare figure [9]). A similar 
comparison of measurements 9 and 21, both taken at a laser power of 25 [iW (right 
panel of figure [8l shows no difference in the ^-profiles within the statistical errors. 




-100 "50 50 100 -100 "50 50 100 

z (pm) z (jjm) 



Figure 8. Comparison of z-profiles of measurement 17 and 18 (left panel) and 9 and 
21 (right panel). Each pair was taken at the same laser power, but different bath 
temperature. For details see text. 



While such a comparison of different measurements seems to be intuitively correct, 
there are several objections against this procedure. First, the measurements we 
compare have quite different values for the spatial temperature and thus different spatial 
extension in the trap. Second, we cannot exclude a priori that, e.g., measurement 18 
has a condensate as well. To overcome these difficulties, we have to look for a way 
to analyse each measurement for itself. Theoretical studies of the thermodynamics of 
an interacting Bose gas of excitons [?] have shown that even in case of a BEC both 
the density of excitons and the spectrally integrated luminescence intensity outside the 
condensate region closely follow a Gaussian distribution reflecting the temperature of 
the exciton gas. Therefore, by fitting a Gaussian only to the wings at larger z of the 
intensity profile, we should have access to the distribution of thermal excitons only. If 
there is no condensate present, this Gaussian should also describe the intensity profile 
in the centre of the trap. On the other hand, any deviation of the measured profile from 
this curve is a clear indication of an additional luminescing component in the trap. If the 
line shape of this component is non-Gaussian, this would clearly indicate the existence 
of a condensate. 

For the measurements shown in figure [7j the results of such a fitting procedure are 
shown in figure [9} The red and blue diamonds are the experimental points, whereby the 
blue diamonds mark those intensity points that have been used in the fitting. The blue 
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Figure 9. Comparison of z-profiles with fits by a Gaussian. Red and blue diamonds: 
total measured profiles, solid blue lines: Gaussian fit of the total profiles, blue 
diamonds: data points used for the fit, green triangles: difference between total profile 
and fit. The numbers in parentheses denote the measurement numbers, cf. figure [6j 
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line gives the shape of the Gaussian. In order to reproduce the obviously different line 
shapes for positive and negative z values due to the Morse-type potential of the trap in 
^-direction, we used for the fit an asymmetric Gaussian function of the form 

S(2) x exp (- ( e( 2o -.>_ + e(l-, Kp) ■ (15) 

For all measurements, the ratio of a + /a_ is obtained as 1.5, in consistence with the shape 
of the potential. The number of data points for the fitting was chosen for each set of data 
such that the average error per point becomes minimum. If we chose less points, the 
statistical error due to noise will increase, if we chose more points, the systematic error 
increases if the profiles contain an additional contribution. The green triangles show 
the difference of the experimental points and the Gaussian fits. In total, the results of 
the procedure substantiate our previous analysis. The measurements at very low power 
(1) and at high bath temperature (20,21), where we expect that there is no condensate 
present, indeed can be fitted completely by a single Gaussian. Measurements 15 and 17 
show clearly a bimodal distribution and thus a strong condensate component. However, 
in measurement 18 a small additional contribution remains. To show systematically 
the results of the fit, we have plotted in figure 10 the ratio of the green component 
to that of the overall intensity. Looking more closely on the results shown in figure 



10] , we can identify at least four sets of data with a strong condensate contribution, 
while most of the other measurements show a contribution below 0.05, which we will 
consider not as significant. The occurrence of a condensate in the power range between 
60 and 300 |J.W is understandable in view of the rather high number of paraexcitons in 
the trap which, according to the rate model (figure [2]) is about 4 ■ 10 7 at 100 |xW. The 
appearance of a condensate at the low power of 5 \xW is rather surprising and requires 
further considerations. Indeed, the theory presented in the next section will show that, 
under the conditions of our experiments (bath temperature 38 mK and no heating of 
the sample), condensation will take place at exciton numbers as low as 4 ■ 10 6 (see figure 



19). 



Here we stress that the luminescence fraction ?7i um shown in figure [TO] is different 
from the fraction of particles in the condensate r\ c = N c /N tot . If f\ um denotes the ratio 
of the luminescence efficiency of excitons in the condensate to that of excitons in the 
thermal cloud (see section 3.2), the relation is given by 

flurcJIc , -, n \ 

Vlum = Z — f • 16 

At present, /i um is neither know experimentally nor theoretically (see section 3.2), 
however, it must be smaller than one because otherwise we would not observe a kink in 



the luminescence dependence on laser power (see figures 11 and 12) so that r] c > r]\ um . 
Our interpretation is further substantiated by the dependence of the totally 



integrated intensity on the excitation laser power which is shown in figure 11 



At low power, the dependence is almost linear, indicating that decay processes by 
exciton-exciton collisions are not important. Actually, the dependence is even slightly 
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1000 



Figure 10. Power dependence of the ratio between integrated intensities of the 
condensate and of the total band. Filled red triangles: measurements 1-3 and 10-17, 
open red triangles: measurements 4-9, blue squares: measurements 18-21. Between 
measurements 3 and 4, and between measurements 9 and 10, there was a change of 
the neutral density filter, cf. figure [6] The shaded band denotes an area where the 
contribution is too small for a significant proof of a condensate. 



1x10- 



1x10 

-t— ' 

| 1x10 ? 

o 



£ 1x10 

CD 
-t— < 

_c 

1x10 E 



1x10" 





10 



100 



1x10 



1.3 



1.1 



0.9 £ 

-t— > 

03 

0.7 ^ 

E 
3 

0.5 



0.3 

3 



P L (UW) 



Figure 11. Dependence of the totally integrated intensity (squares) and spatial 
temperature (triangles) on the laser power. The points with blue colour mark 
measurements 18-21. The dashed line gives a power dependence according to P^J with 
6=1.1 which describes the low power regime quite accurately. 



superlinear which might be due to the effect of increasing homogeneous broadening at 
higher exciton numbers. At power levels of about 60 [iW, we clearly observe a kink 
resulting in a weaker slope at higher power levels. This is just the behaviour predicted 
by the theory in section 3.3.3 The critical power at the kink is in full agreement with 
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the conclusion drawn from figure 10, where we observed the onset of condensation at 
powers above 60 uW. This dependence is abruptly changed for the measurements 18-21 
in agreement with the vanishing of the condensate. At the high exciton numbers (under 
these conditions, we expect N ~ 1 • 10 8 according to figure [2]) we have a substantial 
effect of the two-body decay resulting in a sublinear behaviour, see the data marked by 
blue symbols. 
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Figure 12. Dependence of the totally integrated intensity (squares), spatial 
temperature (triangles), and start and final bath temperature during the measurement 
(black bars) on the laser power for two different sets of data. 



We note that the theory presented in the next section attributes such a kink in the 
power dependence of the total intensity of an exciton gas in a trap to the occurrence of 
condensation. Actually, in all of our measurements using different potential traps and 
under different excitation conditions such a kink shows up, provided that the Helium 



bath temperature is below 400 mK. Exemplarily this is shown in figures 12 and 13 



for three different sets of data. In figure 12 we compare two measurements which 
have the same dependence of the Helium bath temperature but differ in their spatial 
temperatures. While at low powers both curves coincide, the curve corresponding to 
lower spatial temperature shows the kink at lower powers. In figure 13 we display the 
power dependence for the case where both the Helium bath temperature and the spatial 
temperature are much higher and the kink occurs at powers one order of magnitude 
larger than in the measurements of figure 12 This dependence of the onset of BEC on 
both the bath temperature and the spatial temperature will be explained by the theory 
presented in the next section. The data in figure [13] furthermore demonstrate that the 
origin of the kink cannot be the Auger like decay of excitons at high densities, as one 
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might suspect. This deviation from linearity in the power dependence already shows up 
below the kink. 
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Figure 13. Dependence of the totally integrated intensity (squares), spatial 
temperature (triangles) , and start and final bath temperature during the measurement 
(black bars) on the laser power. 



3. Theory of excitonic BEC and their luminescence spectra 

After presenting the experimental results, in this section we will provide the theoretical 
description, attemting to explain the main effects found in the luminescence spectra. 

So far, excitons in potential traps have been described mostly by theoretical models 
of non-interacting bosons. Concepts to include e.g. interaction effects have been worked 
out for atomic condensates. Along this line, the thermodynamics of a one-component 
Bose gas has been investigated extensively, see, e.g., references [?, ?, ?, ?]. First 
applications to excitons exist, too [?]. Recently, distinct signatures of a condensate in the 
decay luminescence spectrum of the non-condensed excitons have been predicted using a 
mean- field formalism in local density approximation [?]. In analogy to generalisations for 
multi-component atomic gases [?, ?, ?, ?, ?] and spinor polaritons [?, ?], a generalisation 
to a multi-component gas of interacting paraexcitons and orthoexcitons has been given 
in reference [?]. 

Here, we focus on the thermodynamics of weakly interacting excitons in a trap and 
assume the whole exciton gas to be in thermodynamic equilibrium. This allows us to 
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apply the theory developed in references [?, ?]. 
3.1. Thermodynamics of trapped excitons 

In order to write down the Hamiltonian of the exciton system, the interaction potential 
between excitons has to be modelled. The description of exciton-exciton interaction is 
a long-standing problem (see, e.g., [?, ?, ?, ?, ?, ?]), and so far no satisfying solution for 
the general case has been obtained. For our calculations, we assume a contact potential 
interaction V(r — r') = Uq5{t — r'), where the interaction strength Uq is given by the s- 
wave scattering length a s via Uq = 4na s /M. M is the exciton mass M = 2.6m e in units 
of the free electron mass m e [?]. This leads to the Hamiltonian in second quantisation 
for the grand canonical ensemble: 



H= J Ah V>t(r,t) [A^- + Kxt(r) - ^ V(r,t) 

+ \jd 3 r U ^(r, t)^(r, t)V(r, t)^(r, t) , (17) 

with the trap potential Kxt( r ) and the chemical potential /i. Decomposing the field 
operators in the usual fashion ^(r, t) = $(r) + ^(r, t), with the condensate wave function 
$(r) = (ip(r,t)) = and the operator of the thermal excitons ijj(r,t), one arrives 

at the Gross-Pitaevskii equation (GPE), 



h 2 v 2 

2M 



n + u 



n(r) + n (r) 



|$(r) + C/ m(r)$*(r) , (18) 



and the equations of motion for the thermal excitons, 



dil>(r,t) 
dt 



2M 



cxt I, 



+ Uq 



r) + m(rj 



with the normal and anomal averages n 1 



-li + 2U n(r)U(r 7 t) 



(19) 

the condensate 



density n c = |$| 



{ip^ip) and m 
and the total density n = n c + n T . 
Since the spatial extension of the potential trap is large compared to the typical 
length scales of the interacting exciton gas (in particular with respect to the thermal 
deBroglie wavelength of the excitons), we can use the local density approximation, thus 
treating the excitons as a locally homogeneous system. The equation of motion (19) is 
solved by a Bogoliubov transformation. In the so-called Hartree-Fock-Bogoliubov-Popov 
(HFBP) limit (m — > 0), the quasiparticle energy reads 



£(k,rW£(k,r)2-(f/ n c (r)); 



(20) 



with 



£(k, r) = Ti 2 k 2 /2M + Kxt(r) - ft + 2U n(r) 



(21) 
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Within these approximations, the non-condensate density n T is given by 



n T (r) 



d 3 k 



E(k,r) 



n B (£(k,r)) 



1 



6 E 2 (k,r 



(22) 



Applying the Thomas-Fermi approximation to the GPE (18), i.e., neglecting the kinetic 
energy term, yields the condensate density as 

1 

Uo 

Evaluating equations (22) and (23), n T and n c have to be determined self-consistently. 



n c (r) 



fi - ^xt(r) - 2lW r (r)] 6(/x - V cxt (r) - 2U n T (r] 



(23) 



3.2. Theory of decay luminescence 

Since the optical wavelength of the emission is much smaller than the trapped exciton 
cloud, we apply a local approximation to the emission spectrum [?, ?] as well, which is 
determined by the excitonic spectral function A(r,k,u), 

I(r,u) oc 2nh\S{k = 0)\ 2 5{hu' - fi)n c {r) 

+ J2 \S(k)\ 2 n B {hu' - fi)A{r,k,hu' - fi) , (24) 

with S(k) representing the exciton-photon coupling. The spectral function is given by 



the quasiparticle spectrum in (|20|): 
A(r, k, uj 



2ixh 



u 2 (k, r)5(hu - E(k, r)) - v 2 (k, r)5(hu + E(k, r)) ,(25) 



where u(k, r) and v(k, r) are the Bogoliubov amplitudes, 

1 



u(k, r) 



v(k, r) 



(£(k,r)/£7(k,r) + l) 
(£(k,r)/E(k,r)-l) 



(26) 



In Cu 2 0, the decay of ground-state paraexcitons in the yellow series is optically 
forbidden. However, due to the applied strain, a mixing with the green series takes 
place and the decay becomes weakly allowed [?]. The paraexcitons decay directly, 
whereby momentum conservation requires that only excitons with the same momentum 
as the emitted photons are involved. This zero-phonon decay can be treated by setting 
uj' = uj — E g x/h with the excitonic band gap E gX and |>S(k)| 2 = Sg5(k — k ) with k 
being the wave vector of the intersection of the photon and exciton dispersions. 

Due to the form of S(k), a condensate of paraexcitons in Cu 2 (energetically in the 
ground state, i.e., k = 0) should not contribute to the luminescence spectrum. However, 
this statement holds rigorously only for homogeneous, infinitely extended systems. Due 
to the finite size of the condensate, it is spread out in k-space, and a weak contribution 
to the luminescence should be expected. Taking the condensate as a classical coherent 
emitter, the radiation follows from classical Maxwell equations as the Fourier transform 
of the polarisation [?] at the wave vector of the emitted photon ko. Since the polarisation 
is proportional to the condensate wave function, the (dimensionless) strength of the 
contribution S c (\S(k = 0)| 2 = S c Sq) should be proportional to the components of the 
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Fourier transform of the ground-state wave function at k = ko- Our calculations show 
that S c is of the order of 10~ 6 to 1CT 7 . Therefore, we will discuss possible effects of 
a BEC in the spectrum considering a weakly luminescing as well as a completely dark 
condensate. 

Furthermore, to account for the finite spectral resolution in experiments, we 
convolve the spectral intensity (24) with a spectral response function of the shape 
exp[— (w/A) 2 ]. Here, A is a measure for the spectral resolution. Using these 
assumptions, the spectrum reads: 

2" 



I(r,u) oc S' c (27r) 3 exp 



huj' — n 



n c (r) (27) 



V A 

+ w 2 (k , r)72 B (£(k , r)) exp [-e 2 _ (a/, k , r)] 

v 2 (k , r)n B (-E(k , r)) exp [-£+(u/, k , r)] , 

with e± (u/, k, r) = (huj' — /x ± E(k, r))/A. As explained in section [2j in the experiment, 
a spectrograph is used to obtain different spectra by integrating over either one or more 
spatial coordinates and/or u. Here, we consider the z-resolved luminescence spectrum 
I(z,u), the ^-profiles of the luminescence I(z), the spatially integrated luminescence 
I(co), and the totally integrated luminescence I tot . The z-resolved luminescence 
spectrum is obtained by imaging a small stripe of width 2 Ay elongated along the in- 
direction onto the entrance slit of the spectrograph, hence integrating over the x- and 
y-direction. The z-profiles are generated by integrating I(z,u) over the energy (to). 
Integrating over all spatial dimensions yields the spatially integrated luminescence I(u>). 
As a fourth option, one can also integrate over r and oj, arriving at the totally integrated 
luminescence J to t, which only depends on the exciton number A" and the temperature 
T. 



3.3. Results 

3.3.1. Luminescence spectrum For the calculations, we used an anharmonic potential 
trap fitted to the experimental of section |2]but with a trap minimum of Vq = 1.04 meV. 
The s-wave scattering length is chosen to be a s = 2.18ae (taken from [?], see also 
[?, ?, ?]) with the excitonic Bohr radius a-Q = 0.7 nm. 

First, we revisit the flat bottom shape of the spectrum discussed in reference [?]. 



In figure 14 we plot the luminescence spectrum for a constant temperature T = 0.6 K 
and three different exciton numbers (left column: = 2.5 ■ 10 7 , middle column: 
= 2.5 ■ 10 8 , right column: = 8.0 ■ 10 8 ). The condensate is expected to remain 



completely dark. The left image of figure 14 shows a thermal spectrum exhibiting the 



typical nearly parabolic shape that was also found in the experimental results presented 
in figures [3] and [7j Increasing the exciton number by a factor of 10 while keeping 
the temperature constant results in the onset of a BEC with a condensate fraction of 
N c /N = 0.05. The shape of the spectrum is altered and develops a flat bottom at the 
chemical potential as reported in reference [?]. Further increasing the exciton number 
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Figure 14. z-resolved luminescence spectrum I(z,u>) at T = 0.6 K with dark 
condensate. The exciton numbers are N — 2.5 ■ 10 7 (left), N = 2.5 • 10 s (middle) and 
N = 8.0 • 10 8 (right). Condensate fractions are N c /N = (left), iV c /iV = 0.05 (middle) 
and N c /N = 0.50 (right). Trap minimum is at Vq = 1.04 meV and e — Tllu — E g x — Vq. 



leads to a more pronounced flat bottom and an energetic shift of the spectrum with the 
chemical potential. Additionally, a new contribution to the thermal spectrum below the 
chemical potential, arising from the t> 2 -term in equation (25), becomes faintly visible. 
These effects are linked to the occurrence of a BEC and would indicate its existence 
even without direct emission from the condensate. 

Comparing these predictions with the experimental findings, e.g., figures [3] and [7j 
one has to conclude that the signatures were not observed in the experiment. However, if 
the condensate exhibits a very weak luminescence, the changes in the spectrum predicted 



from theory are much more subtle. This can be seen in the left column of figure 15 



(S c = lCT 6 /3), cf. again figures [3] and [7j In contrast to the case of a dark condensate, 
there is no drastic qualitative change in the spectra from top to bottom. Also, the shift 
on the energy axis appears to be smaller and could well be an interaction effect of thermal 
excitons. Therefore, the contribution from the condensate hides the flat bottom as well 
as the shift on the energy axis, without being a clearly visible delta-shaped peak, as 
expected for a fully contributing condensate. This explains why a conclusive detection of 
a condensate via spectral signatures requires a very careful analysis of the experiments. 

The spectra in the left column of figure [15] consist of different contributions. To 
illustrate this, we plot the results for the spatially integrated luminescence I(u) in the 
middle column of figure 15, cf. the measurements shown in figures [3] and [5| The upper 
panel shows a spectrum in the non-condensed case at low densities. Its shape is given 
by a Bose distribution convolved with the spectral resolution of the spectrometer. The 
middle panel shows a spectrum at higher densities, where interaction effects are already 
important and a very small condensate contribution already occurs, nearly invisible 
in the total curve. Compared to the previous case, the peak becomes narrower and 
shifts to higher energies, while the tail remains qualitatively the same. The bottom 
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Figure 15. z-resolved luminescence spectrum, z-integrated luminescence, and z- 
profiles of the luminescence with weakly luminescing condensate (S c = 10 _6 /3) for 



the same parameters as in figure 14 Left column: luminescence spectrum I(z,uj). 
Middle column: z-intcgrated luminescence I{<jj). Right column: z-profiles I(z). Red 
curves (full and dashed) denote the contributions from thermal excitons, the blue dash- 
dotted curve gives the condensate contribution, and the total emission is shown as a 
dashed black line. 



graph shows the case with a condensate fraction of 0.5 including a distinct condensate 
contribution (dashed-dotted blue) which has a Gaussian shape and is centred at the 
chemical potential. The contribution from the thermal excitons consists of the u 2 part 
of the spectral function (solid red) and v 2 part of the spectral function (dashed red), 
compare equation (27). The black line represents the sum over all contributions. In 
contrast to the very weakly condensed case (middle), the contribution related to the 
v 2 term appears below the chemical potential (dashed red). Except for this additional 
contribution, the total emission is only slightly shifted towards higher energies and has 
a higher maximum compared to the middle graph in figure 15 The width of all the 



Gaussian like peaks in figure 15 is directly given by the spectral resolution A = 18 [ieV. 



3.3.2. Spatially resolved luminescence The results for the spatially resolved 
luminescence I(z) are presented in the right column of figure 15 and should be compared 
to the experimental results shown in figures [3j [3], |HJ and [9j 

Without any contribution from the condensate, one would expect the solid red 
line. The non-condensed case follows a Gaussian shape (top graph), while the onset 
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of the BEC leads to a deformation in the form of a plateau (middle and bottom 
graph). However, taking a weakly luminescing condensate (dashed blue) into account, 
the total emission (black dashed) looks Gaussian like again, masking the signature of 
the condensate (middle and bottom graph). 




Figure 16. Totally integrated luminescence /tot for T = 0.5 K and different fractions 
of condensate contribution. Change of -/V cr it (position of kink) with temperature shown 
in the inset. 



3.3.3. Totally integrated luminescence Integrating over all variables (r and oS) yields 
the totally integrated luminescence 7 to t, which had been considered in the experiments 
in figures TTJ 12 , and 13 In figure 16 , we show the numerical results for 

-^tot & 

function 

of the exciton number iV for a fixed temperature of T = 0.5 K. Remember that N is 
connected to the excitation laser power P L , cf. the results of the rate model, figure 
[2j For small N, the totally integrated luminescence increases linearly with the exciton 
number until a critical value is reached. At the critical exciton number, the curve has a 
kink and continues with a weaker slope afterwards. The additional contribution from a 
weakly luminescing condensate does not alter the behaviour qualitatively in this case, 



as it can be seen from the other curves in figure 16 As shown in the inset of figure [16} 
the critical exciton number shifts with the temperature approximately as expected with 
T 3 , equation (Jl]). Obviously, here the interaction has no drastic effect on the behaviour 
known from an ideal gas. 

However, in the experimental situation, it is merely impossible to measure over 
orders of magnitude of the particle number without changing the exciton temperature. 
The temperature will rise as a result of the energy introduced into the system by 
the laser. Therefore, experimental points may lie anywhere in the N-T plane. The 
theoretical results for the totally integrated luminescence as a function of exciton number 
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and temperature are shown in figure 17 Depending on how fast the temperature 
increases with exciton number, it is very well possible to never cross the phase boundary 
and observe the kink in the totally integrated luminescence. 

3.4- Discussion 

In the previous section, we presented different spectral signatures of a non-emitting 
condensate in the thermal emission of the excitons in global thermal equilibrium. The 
five main signatures are: (i) the formation of a flat bottom in the luminescence spectrum 
I(z, u), (ii) the shift on the energy axis of the spectrum I(z, u>), (iii) a deviation of the z- 
profiles I(z) from the Gaussian shape, (iv) the appearance of the v 2 term in the spatially 
integrated luminescence (v) the kink in the totally integrated luminescence J tot at 

iV crit . 

If the contribution from the condensate is much weaker than what we estimated 
(i.e., S c of the order of 10 _6 ...10~ 7 ), signature (i) should be experimentally visible. This 
is neither the case in the current experiment nor have other experiments seen this effect. 
If the contribution is much stronger though, one should expect the delta-shaped peak 
as predicted by the standard theory [?, ?, ?]. This substantiates our estimation for S c . 

However, if the condensate contribution is comparable to that of the thermal 



excitons, the flat bottom in (i) can be masked as shown in figure 15 The changes 
in the shape are very subtle and probably not detectable in an actual experiment. The 
shift on the energy axis (ii) is also altered by the condensate contribution. Additionally, 
there are numerous other effects that can change the energetic position of the spectrum, 
e.g., a background plasma of electrons and holes [?, ?]. Therefore, this shift might 
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not be a good indicator for the onset of the BEC. The deformation of the ^-profiles (iii) 
would also be masked by the condensate contribution. Here, the latter one complements 
the thermal luminescence to just again a Gaussian profile. Contrarily, in the measured 
profiles (figure [9]), the condensate adds a contribution to the thermal Gaussian profile. 
We will revisit this point later. 

Without a condensate, the spectra can be fitted with a renormalised Bose 
distribution convolved with the spectral resolution of the spectrometer. With a 
condensate, the v 2 peak (iv) causes a characteristic deformation of the low energy flank 
of the spatially integrated luminescence. For a dark condensate, however, it would 
be a free standing peak separated from the rest of the spectrum. Both cases should be 
detectable in an experiment though. However, this signature (iv) as well as the energetic 
shift (ii) are very subtle and are probably masked by noise, compare figures [3] and [7j 

The most promising signature of condensation would be the kink in the totally 
integrated luminescence (v), as it is not altered qualitatively by the condensate 
contribution. In this case it only changes the slope after the kink, but does not hide 



the kink itself. Although in the figures 11 T3mo isotherms are plotted - the exciton 



temperature rises with excitation power - there is obviously a kink in the measured 



integrated intensity. Looking at figure 17 a path on the surface will exhibit a kink even 
with increasing temperature, if the ridge is crossed. Therefore, we can relate the kink 
in the experimental figures, at least qualitatively, to the occurrence of a condensate. 

3.5. Excitons in local equilibrium 

As stated above, the theory derived in section [3] assumes that the exciton gas in the 
trap is in global thermodynamic equilibrium. The experimental results presented in 
the previous section show, however, that this is obviously not the case and that the 
exciton temperature is the key quantity to detect deviations from global equilibrium. 
The conclusions are that (i) the spectral (or spatial) temperature is not equal to the 
exciton temperature and (ii) the assumption of global equilibrium, therefore, must be 



wrong, cf. the discussion in section 2.3.1 



In contrast to the latter statement, the luminescence spectra calculated by the 
equilibrium theory reproduce the measured spectra qualitatively quite well, cf . figures [3] 



and 15 Thus it seems to be reasonable to abandon the global equilibrium assumption, 
but to keep local equilibrium. In this case, the thermally excited excitons have still an 
equilibrium (Bose) distribution, but with spatially varying temperature and chemical 
potential. The space dependence of the latter quantities is in principle unknown, but 
from the experiment we can conclude that the temperature should not vary over the 



dimension of the trap, see section 2.3.1 The only criterion for the spatial dependence 
/Lt(r) is that the spectra must be reproduced. Therefore, we demand that the spatially 
integrated spectrum I{u) (at least its high-energy tail) follows a Bose distribution with 
the spatial temperature and some formal constant chemical potential jx which is used 
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to fix the particle number. This leads to 



1 - 



T 

-*- R 



(V eyLt (r) + E + 2U n T (r)) + ^fi 



(28) 



which allows to directly use the theory of section 3.1 also in this nonequilibrium situation. 




Figure 18. Spatially resolved luminescence spectrum (left) and ^-profiles of the 
intensity (right) for N = 10 8 and T s = 0.6 K. Upper row: exciton temperature 
Tx = 0.25 K (uncondensed case), lower row: Tx = 0.15 K (condensed case, r\ c = 0.2). 
Solid red lines: total profiles (thermal + condensate), solid blue: only thermal 
contribution, dash-dotted blue: Gaussian fit of the total profiles, blue crosses: data 
points used for the fit, solid green: difference between total profile and fit. The ratio 
between integrated intensities of the condensate and of the total band is ?7i um = 0.07 



In what follows, we apply the local equilibrium theory to two typical experimental 
situations. Thereby we fix the total particle number N = 10 8 and the spatial 
temperature T s = 0.6 K and vary the exciton temperature. The condensate contributes 
to the luminescence with S c — 4 • 10~ 7 . The spatially resolved spectra for both cases are 
displayed in figure 18, left column. As expected, in the uncondensed case (T x = 0.25 
K, upper row), the spectrum just follows the external potential. The condensed case 
(Tx = 0.15 K, lower row) with a condensate fraction r] c = 0.2 looks qualitatively not very 
different. In particular, due to the spatial variation of the chemical potential according 
to equation (28), the flat bottom of the spectrum (cf. figure 14) disappears. 

In the condensed case (Tx = 0.15 K), the z-profile of the luminescence exhibits a 
clearly non-Gaussian shape (figure |i~8| right column, lower panel, red line). For the not 
much higher temperature of T x = 0.25 K (upper panel), the shape is approximately 
Gaussian. Therefore, we fit the thermal component of the z-profile for the condensed 
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case by an asymmetric (due to the potential asymmetry) Gaussian. In order to 
exclude the condensate contribution, we omit in the fitting procedure the data points 
in the centre of the trap where the condensate is situated. The size of the excluded 
area is determined by minimising the error of the fit. The fitting result is given by 
the dash-dotted blue line in the lower right figure. It basically follows the thermal 
profile component. The small deviations from the Gaussian shape are caused by 
the potential anharmonicity and by the renormalisation of the potential due to the 
interparticle interaction. Thus, the procedure to extract the thermal component from 
the experimental ^-profiles shown in figure |9j where a Gaussian fit has been applied, 
too, seems to be justified. The difference between the total ^-profiles and the Gaussian 
fits is given by the green lines. It basically reflects the condensate contribution. For 
the condensed case we obtain a luminescence fraction r}\ ura = 0.07 which would set 
/lum = 0.25. 

The two cases depicted in figure 18 can be compared to the measured z-profiles 
shown in figure [9j The "uncondensed case" (upper row in figure 18) obviously 
corresponds, e.g., to measurements 1 and 21, where the difference between data and 
fit is just noise. On the other hand, the "condensed case" (lower row in figure 18 ) finds 
its counterparts, e.g., in measurements 3 and 17. In each striking qualitative 

agreement is found. This corroborates our explanation of the experimental findings by 
the occurrence of a excitonic BEC substantially. 

Figure 19 shows the critical particle number for Bose-Einstein condensation in 
dependence on exciton temperature Tx and spatial temperature T s . The global 
equilibrium case is represented by the diagonal Tx = T B . If the excitons are excited 
outside the trap, their spatial profile will be at first always broader than in equilibrium. 
Therefore, T s < Tx can be ruled out in the current experiment. Obviously, for a given 
exciton temperature, the critical number is smallest for global thermal equilibrium and 
increases with increasing spectral temperature. Even at low powers, i.e., for quite low 
exciton numbers, a condensate is possible if the exciton temperature is close to the He- 
bath temperature, cf. the situation shown in figure [9j measurement 3, where the spatial 
temperature is 0.35 K and the exciton number is 4- 10 6 which just corresponds to the 
critical number at Tx = 0.04 K. 

With the help of this figure, we can systemise the experimentally observed power 
dependence of the totally integrated intensity shown in figures 11-13 This is done 
by looking at which exciton temperature Tx the critical number (corresponding to the 
critical power and the respective spatial temperature) crosses the critical surface. These 
points are plotted as coloured rectangles in figure 19 The values for T x give a very 
systematic picture. With increasing critical power, both the spatial and the exciton 
temperature rise. Thereby, Tx is always larger than the bath temperature, which is 



consistent with the expected crystal heating, compare Appendix C 



Figure [19] also compares the data for the measurements 17 (blue star) and 18 (red 
star) of figures [7] and [9] showing that the conditions in measurement 17 are such that we 
are well above iV cr i t while for measurement 18 we are just touching the critical surface. 
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4. Conclusions and Outlook 

We have presented a series of experiments that investigates the luminescence from 
paraexcitons in cuprous oxide after cw-excitation. The excitons have been confined in 
a potential trap with the host crystal immersed in liquid helium at temperatures down 
to 38 niK. The spatially resolved luminescence spectra do not show notable variations 
when varying the laser excitation power from below 1 (xW to above 1 mW. However, we 
observed a number of characteristic changes in the spatial profiles of the luminescence 
and in the totally integrated intensity: In the spatial profiles, a new component in the 
centre of the trap occurs at intermediate laser powers which cannot be fitted by a thermal 
Gaussian distribution. This component vanishes at very low power and at very high 
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Figure 19. Critical particle number N clit for Bose-Einstein condensation vs. exciton 
temperature Tx and spatial temperature T s . The diagonal Tx = T s corresponds to 
global equilibrium; Tx > T s is, for excitation outside the trap, physically meaningless. 
The blue coloured rectangle corresponds to the critical power of 30 (J.W for the data 



shown in figure 12 (blue symbols), the red coloured rectangle corresponds to the critical 



power of 60 u.W for the data shown in figure 12 (red symbols), the black coloured 



rectangle corresponds to the critical power of 600 |J.W for the data shown in figure [13} 
and the green coloured rectangle corresponds to the critical power of 78 (J.W (data not 
shown) . The magenta coloured dot corresponds to the critical power of 40 p.W from the 
measurements shown in figure [XT] The stars denote the conditions of measurements 17 
(blue) and 18 (red), the black diamond denotes the critical number for measurement 
3, all of figure [7J 
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power levels. At the mentioned intermediate power levels, we see deviations from the 
linear dependence of the totally integrated intensity on the laser power: Most notably a 
well defined kink appears in the slope of the curves. From the linear dependence of the 
total intensity on laser power, we conclude that two-body decay processes of the excitons 
do not play any significant role. This is in agreement with previous measurements under 
pulsed excitation in the same samples [?]. 

With increasing laser power, we observe a slight increase in the effective temperature 
that is inferred from the high-energy tail of the luminescence spectra from 0.35 K up to 
a maximum of 1 K. This spectral temperature does only reflect the spatial distribution 
of excitons in the trap and is not identical to the local temperature which determines 
the energy distribution of the excitons at each point in the trap. The local temperature 
will be determined by exciton relaxation processes and is expected to go down almost 
to the bath temperature at low excitation power. 

Recently, measurements [?] reported a strong heating of the exciton gas at high 
pump powers and claimed that this effect originates from a relaxation explosion of 
excitons when a transition into a Bose-Einstein condensate takes place. We observe a 
similar heating under high excitation powers but can definitely rule out the existence of 
a BEC in this high power range. 

In order to understand the observed features of the luminescence spectra, we 
theoretically analysed the thermodynamics and the luminescence properties of excitons 
in a potential trap. Thereby the excitons are described as an interacting Bose gas 
in the framework of a Hartree-Fock-Bogoliubov-Popov approach. Already under the 
assumption of global thermal equilibrium, this theoretical approach suggests that the 
observed kink in the totally integrated intensity signals the transition into a Bose- 
Einstein condensate of trapped paraexcitons. Taking into account the specific non- 
equilibrium situation in the trap, the theory also consistently describes the characteristic 
changes in the spatial profiles of the luminescence: While the theory predicts almost no 
changes in the spatially resolved spectra, it allows to identify the additional component 
in the spatial profiles as due to a weakly luminescing condensate. 

To conclude, we have presented strong evidence, both from experiment and theory, 
that at ultracold temperatures in the range of 100 mK, paraexcitons in Cu 2 undergo 
a transition into a Bose-Einstein condensate. 

Nevertheless, further experimental investigations are necessary to prove the 
existence of a BEC of excitons in cuprous oxide beyond any doubt. For example, 
direct measurements of the lattice temperature and of the local exciton temperature 
will provide a better understanding of the thermodynamics. Thereby, it might be 
intriguing to apply spatially resolved Brillouin scattering since the energies of the 
phonons involved are comparable to the thermal energy. The local exciton temperature 
would be detectable via infrared absorption of the IS — 2P transition [?, ?]. The 
most important aspect, however, must be the direct proof of the macroscopic coherence 
of the condensate by interferometric methods or by means of intensity correlation 
measurements, both being under way. 
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From the side of the theory, one can try to improve the exciton density calculations 
by using more advanced approximations, e.g., by solving the Gross-Pitaevskii equation 
exactly instead of using the Thomas-Fermi approximation. It also seems to be necessary 
to improve the luminescence theory in order to account for the inhomogeneity of the 
system more rigorously, e.g., following the ideas of reference [?]. Moreover, one should 
also include the spectral broadening due to exciton-exciton interaction. This would 
require the inclusion of higher order correlations beyond the Hartree-Fock-Bogoliubov- 
Popov approximation by calculating the densities and the spectral function on the level 
of the Beliaev approximation [?]. Recently, we have shown that, already for the thermal 
excitons, a more realistic description requires to take the lifetime broadening of the 
exciton states into account [?], which relaxes the strict wave vector conservation in (24). 

Quite recent results by Naka et al. [?] show that the Auger-decay produces a 
significant number of free electrons and holes which are also captured by the trap (see 
figure [2]). This should lead to a shift and a broadening of the exciton states by the 
surrounding electron-hole plasma [?, ?, ?, ?]. An inclusion of this "plasma damping" 
requires to include collisions of the excitons with charged fermions in the theory. This 
is clearly beyond the scope of the theory of weakly interacting bosons. 
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Appendix A. Strain Hamiltonian for electron-hole states 



For the derivation of the energy shifts of electron-hole pairs we consider that the top Ty" 
valence band states can be written as [?] 

*+l/2 = ' 



/2 



1 

7t 



(KT 2 - K?K + 2Y 2 ~% 



(KT 2 - Y 2 2 )(3 V - 2Y 2 1 a v 



(A.l) 
(A.2) 



Here Y 2 m denote the spherical harmonics and a CjV , /3 C)V the electron spin functions of 
valence and conduction band. By inspection of table VI of reference [?] we immediately 
see that the paraexciton state 



$12 = $ ys^| (Vi/ 2 ac + *+i/ 2 A 



(A.3) 
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with $ys being the envelope function of the yellow IS state. We immediately see that 
the paraexciton state has the same behaviour as the r^" band states. In all matrix 
elements of the strain Hamiltonian, all terms containing the electron-hole exchange are 
missing. Since the interaction with the Tg states does not change, the electron-hole 
pairs show the same behaviour under strain, i.e., they perceive an effective trapping 
potential similar to the paraexcitons. 

Appendix B. Exciton relaxation 

A central problem in the dynamics of excitons at milli-Kelvin temperatures is the 
relaxation and thermalisation by a contact to a bath of thermal phonons, as the 
acoustical modes will freeze out [?]. Therefore, we have simulated the relaxation of hot, 
laser excited excitons by assuming interactions with longitudinal acoustic phonons and 
Auger-like two-particle decay in a potential trap by integrating the Boltzmann equation 
[?]. Here we present additional results for the homogeneous situation including elastic 
exciton-exciton scattering, where the model follows closely that described in [?, ?, ?]. 
The resulting system of differential equations was integrated using as initial distribution 
a Gaussian of width 0.1 meV centred at = 5 meV. The initial exciton density was 
assumed to be no = 1 ■ 10 15 cm -3 . For the elastic exciton-exciton scattering cross section 



a = 50 nm 2 [?] was taken. The main results are shown in figure Bl, where the effective 
temperature of the exciton gas obtained by fitting the distribution function by a Bose 
distribution (see inset) is plotted for different lattice temperatures as a function of time. 



Appendix C. Thermal behaviour of the dilution refrigerator and the sample 

The 3 He/ 4 He dilution cryostat uses a mixture of 3 He and 4 He for the cooling process. The 
coldest region is inside the mixing chamber, consisting of a 3 He-rich (100% 3 He) and a 
3 He-poor (6.4% 3 He) phase separated by a phase boundary, the basic thermodynamics 
of which is well known [?]. In the following analysis, we assume that the sample is 
immersed inside the dilute 3 He/ 4 He mixture and that the resistor measuring the bath 
temperature is placed in between the sample and the phase boundary in the mixing 
chamber. 

To derive a connection between sample temperature and measured bath 
temperature we have to consider the whole system as a thermal network as shown 



in figure CI The incoming laser beam with the power Pl hits the sample which is 
at a temperature Ts and has a heat capacity C v (Ts)- It is surrounded by the dilute 
3 He/ 4 He mixture (temperature T# e ), which itself is connected to the mixing chamber at 
a temperature Tm with heat capacities C v (Tu e ) and C v (Tm), respectively. The mixing 
chamber is cooled with a power P C (T M ) depending on its temperature. The heat 
conduction between the different compartments is represented by the heat conductivities 
Ai and A2. Applying the continuity equation for the energy flow, the network is described 
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Figure Bl. Cooling behaviour of laser excited paraexcitons for different bath 
temperatures. Plotted are the effective temperatures of the distribution function 
obtained by fitting with a Bose distribution. The cooling curve for a lattice temperature 
of 0.05 K is identical to that of 0.1 K. The inset shows as an example the distribution 
for a lattice temperature of 0.3 K, an exciton density of n = 1 • 10 15 cm~ 3 at t = 100 
ns from which a temperature of T e g — 0.33 K is derived. 
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Figure CI. Thermal network describing sample and mixing chamber in the 3 He/ 4 He 
dilution cryostat. For explanation, see text. 



by the following system of equations: 

jrri 

a(T s ) ■ ~ = F heat (P L ) - Ai ■ (T s - T He ) , (C.l) 

CV(Thc) • = Ax • (T s - T Hc ) - A 2 ■ (T Hc - T M ) , (C.2) 

C7 V (T M )~ =X 2 -(T Ue -T u )-P c (T M ). (C.3) 

Here Fh ea t gives the connection between the laser input power and the heat generated 
in the sample due to the excitonic relaxation and decay processes (compare section 



2.2). For the parameters given in table 1, the results are shown in figure C2 We can 
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approximate the dependence by a function of the form 

Fhcat = aP L + f3P XX (yjl + P L /Pxx - l) (C.4) 
with the fit parameters a = 0.045, ft = 0.007, and Pxx = 16.67 (J.W. 




according to equation |C.4| 



Heat capacities 

Regarding the heat capacity of the sample, it is necessary to take into consideration 
that not only the crystal is heated by the laser, but, due to the good thermal contact 
supplied by the pressure on the sample, also the sample holder, which consisted of pure 
titanium. Therefore, we assumed the temperatures of sample and sample holder to be 
the same. The variation of the heat capacity of the sample CU2O with temperature 
is well approximated by the Debye model as a T 3 temperature dependence. A fit of 
measured data [?] by 

Cv Cu2 o=A-T B , (C.5) 

results in the parameters A = 8.06 ± 0.05 J/(mol K 4 ) and B = 3.0 ± 0.05. For 
the specimen holder, we consider only the superconducting state below the critical 
temperature (T CT . = 0.4 K [?]), restricting the following analysis to the interesting case 
Ts < 0.4K. Here, the heat capacity is composed of an electron and a phonon part, 
C VTi = C Vc + C Vph . While the phonon part is negligible, the electron part is given by 

CV Ti = C Vc = C Ti 7Ti T c exp ^ 1 ^ Tc ^j , (C.6) 
with C*Ti = 9.17 and the Sommerfeld constant 7ti = 3.3mJ/(mol- K 2 ) [?]. 
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Experimental data for the heat capacity of the helium mixture [?] for different 
concentrations of 3 He in 4 He were interpolated to a concentration of X3n e = 6.4% and 
fitted resulting in a dependence 

a(T He ,M) = D ■ (^) E + F, (C.7) 
with D = 0.417 J/(molK), E = 6.117 and F = 0.7 J/(molK). 

Cooling power 

For the description of the cooling power Pq of the dilution process, where the 3 He- 
atoms flow from the 3 He-rich (concentrated) in the -poor (dilute) phase, we consider the 
enthalpy of the system. By assuming that pure 3 He flows through the heat exchangers 
into the mixing chamber with the dilute phase, the cooling power is given in accordance 
with [?] by: 

P c = n 3 -[H d (T)-H c (T)} (C.8) 
= ?V[95-T M -11-T^], (C.9) 

with the molar flow rate n 3 , the temperature in the mixing chamber Tm and the 
temperature behind the heat exchanger TV. To obtain the unknown parameters of 



equation (C.9), we take as calibration points the lowest temperature reached without 
any heat load Tm = 20 mK and Pq = 100 |aW at Tvi = 100 mK. From these data we 
obtain n 3 = 1.25 ■ 10~ 4 mol/s and Tw = 0.115 K. From this cooling power we have 
to subtract the heat load Pq due to the cryostat windows, which lead to a minimum 
temperature of Tm = 38 mK corresponding to a power of Pq = 13.4 \xW. 

Determination of the temperatures in the cryostat system 

To obtain the temperatures of the different parts of the thermal network from the system 
of coupled differential equations, it is necessary to specify the amount of substances in 
the system. 

The CU2O specimen is a cube with an edge length of 3 mm, a molar mass of 
M = 143.09 g/mol and a density of p = 6.48 g/cm 3 . Hence, the amount of substance 
inside the cube is n Cu20 = 1.398 • 10~ 3 mol. From the measured weight of the titanium 
sample holder of m = 26 g and a molar mass of M = 47.9 g/ mol, we obtain = 0.5 mol. 
The amount of substance in the helium mixture can only be estimated. In our system, 
a 15 % helium 3 He/ 4 He mixture of approximately V = 400 1 circulates under a pressure 
of p = 0.8 bar. That leads to an amount of helium mixture of n3He/*He — 14.3 mol 
from which only about two thirds are in the mixing chamber. For our calculations, we 
assume that 1 mol is in the bath and 2.75 mol is in the mixing chamber. Then the only 
quantities which have still to be specified are the heat conductivities Ai 2- This was done 
by adjusting the solutions of the coupled differential equations to different measurements 
until all data could be described by the same set of parameters. As examples we show 
in Figure C3 the results for the series of measurements from Figure [5] and [6} 
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Figure C3. The numerical solution of the system of differential equations expresses 
the variation of the different temperatures (cf. equations (C.2 |-([C3|) with the 
measurement time for Ax = 6-10~ 4 W/K and A 2 = 27-10~ 4 W/K. The solid red 
framed bars represent the temperature of the sample Tg, the blue dotted curve 
shows the temperature of the surrounding helium bath Tjj c and the blue dashed line 
describes the temperature in the mixing chamber Tyi. The solid blue framed bars are 
the representation of the measured incoming laser power P^. The measured helium 
temperatures are shown as black diamonds. 



The calculations were done with Ai = 6- 10" 4 W/K and A 2 = 27- 10~ 4 W/K and 
show an almost quantitative agreement with the experimental data, despite very large 
changes in laser power (solid blue framed bars) during the measurements, which indicates 
the correctness of our model. 
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